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Abstract 

We study a simple system of two hyperbolic semi-linear equations, 
inspired by the Einstein equations. The system, which was introduced 
in pQ, is a model for singularity formation inside black holes. We show 
for a particular case of the equations that the system demonstrates a 
finite time blowup. The singularity that is formed is a null singularity. 
Then we show that in this particular case the singularity has features 
that are analogous to known features of models of black- hole interiors 
— which describe the inner-horizon instability. Our simple system may 
provide insight into the formation of null singularities inside spinning 
or charged black holes. 
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1 Introduction 

This research examines the singularity formation in a system of equations 
motivated by Einstein's equations. This study may lead to a better under- 
standing of the internal structure of black holes and more specifically, the 
singularity created within them. The system of equations studied in this 
work is motivated by the former studies of the Reissner-Nordstrom (RN) 
based models [T] . It is a continuation of the work in p] where the system of 
equations and their physical motivation were presented. Here we present a 
new approach to the study of the model and prove the conjecture presented 
in [1] for a simplified form of the equations. 

The mass inflation model presented in [3] demonstrates the genericity 
of null singularities inside realistic black holes. In [1] it was conjectured 
that a simpler dynamical system - an "active ingredient" of the Einstein 
equations is capable of producing black hole-like configurations of generic 
null weak singularities subjected to certain characteristic features. It was 
also shown that the evolution equations in Reissner-Nordstrom-de Sitter 
(RNDS) spacetime with two intersecting null fluids can be written in the 
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form: 



F(R) is defined as: 
where 



R )UV = e s F(R), (1) 
s, uv = e s F'(R). (2) 

F(R) = -h'{R) (3) 



R = r 2 , e s = rf, (4) 
for the following metric in double- null coordinates (u,v): 

ds 2 = -2f(u,v)dudv + r 2 (u,v)dQ 2 , (5) 

and 



,1 , 2Q 2A 
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h(R) = 2R2 -Am + — 1 R*, 

R2 3 

where m,Q and A are the mass, charge and the cosmological constant in 
the RNDS solution. The conjecture was formulated in pQ for a general 
function h(R) which has three roots h(Ri) = (Fig. [T]). The roots of h(R) 
correspond to the locations of the "horizons" in a flux-free solution (defined 
below). R = Ri, R 2 , R3 correspond in the RNDS model to the cosmological, 
event and inner horizons respectively. The origin of axes is placed at R = R a . 
In addition we define 

ki = \h'(Ri)\ (i = 1,2,3), 

as a generalization of the surface gravity. In this work we study the sin- 
gularity formation for a simple type of h{R) - a saw-tooth function. We 
will show that the formed singularity and its properties depend only on the 
surface gravities at the horizons and their relative values. 

Our semi-linear system is accompanied by two conserved fluxes: 

e s (e- s R, u ), u = *(«), (6) 
e"{e-'R, v \ v = *(v), (7) 

where \l/(it) and $>(v) are arbitrary functions which stand for the pertur- 
bations of the flux- free solution ^f(u) = $(y) = 0. These functions are 
proportional to the energy density in the outflux and influx respectively. 
Eqs. (|6|7jl can be derived from our semi-linear system (see [1]). Too strong 
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Figure 1: The function h(R) - this function can be any function which has three 
roots: R\ , R2 and R3 . R a is the point where we will place the origin of the axes in 
the next sections. 

null fluids may ruin the singularity structure. Thus, we consider "small" 
ty(u) and $(v). Later we will define exactly what we mean by "small". 

In addition there exists a gauge freedom (u — > u(u),v — > v(v)) for which 
s(u,v) is changed according to the following rule: 



There are two quantities which play an important role in the conjecture and 
are related to the mass function and the null fluid fluxes: 

• e~ s R tU i? „ - A gauge invariant quantity. 

• e~ s R jU (res. e~ s -R „) - Invariant under gauge transformations in u 



The paper is organized as follows. In section 2 we study the singularity 
formation in the flux- free case using the method of characteristics. This 
method will allow us to study the singularity formation in the non-linear 
case (with two null fluids - perturbations of the flux-free case) as well. For 
this purpose we introduce in section 3 the general form of characteristic 
equations. In the discussion (section 4) we review the conjecture about the 
singularity formation (which was presented in [I] ) and indicate the particular 
case which will be studied in the rest of the paper - h{R) as a saw-tooth 
function. In sections 5 and 6 we prove the conjecture for a saw-tooth function 
under few reasonable assumptions. 

2 Singularity in the flux-free case and the method 
of characteristics 

Our semi-linear system (|H2p for the functions R(u, v) and s(u, v) is a system 
of second order hyperbolic PDE. To obtain a unique solution inside the 





(res. v). 
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positive quadrate (u > 0, v > 0) we should have four known functions 
on the boundary of the region: R(u, 0), R(0,v), s(0,v) and s(u, 0). This 
type of initial value problem is usually called a Goursat Problem(see for 
example [5]). Notice that the initial data is given on the characteristic 
curves of the equations. The Goursat problem for linear hyperbolic PDE 
is a well-posed problem and has a global unique solution. Our system of 
equations is semi-linear, so we have only local existence and uniqueness in 
general. This local existence and uniqueness can be proved using successive 
approximations as in the linear PDE case. 

In addition to the initial data we have a gauge freedom in the coordinates 
as our set of equations is invariant under gauge transformations. We can 
choose any function for and || (recall the gauge transformation ([8])) . 
This freedom enables us to set a gauge in which s(0, v) = and s(u, 0) = 0. 
We will refer to this gauge later as the "standard gauge" . In this gauge we 
still have to set only two functions: R(u,0) and R(0,v). Relating to these 
two functions, we divide the initial value problem into two cases: 

1. The linear case - R(u,0) and R(0,v) are linear functions. This is the 
case when we do not have perturbations (flux- free). 

2. The non-linear case - The two functions R(u, 0) and R(0, v) are non- 
linear. 

In the standard gauge, for the flux-free case, eqs. (|6l7p on the boundary 
v = and u = are 

R,uu = (v = 0), (9) 
R m = (u = 0). (10) 

Therefore, on the u axis R is linear with respect to u and on the v axis R 
is linear with respect to u. It is consistent with eqs. (|9ll0p to assume that 
R(u, v) is monotonous on the axes: 

f R(u, 0) = R a - bu / , n\ -in 
\R(0,v)=R a + av, ^ >0 ) W 

where R a as noted in Fig. [1] is placed at the origin of axes. Equating the 
expression e~ s R jV R tU (which is gauge invariant) in the static gauge to the 
same expression in the standard gauge at the point R = R a , we obtain: 

- h(R a ) = -e~ s R 2 x = e~ s R )V R )U = -ab. (12) 
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Eqs. (|6|7p in the flux- free case takes the following form: 



e s (e- s R v ), 



0, 
0. 



(13) 
(14) 



Our goal is to derive a first order PDE whose characteristic curves will 
determine the behavior of the solution of the original semi-linea system. We 
start by integrating eqs. (j!3|14p : 



R,u = e s (j){v) 
R iV = e s ip(u) 

(y),ijj(u) are defined below. After differentiation, we have: 



R,uv 
R.vu 



e s 4>'(v) + e s s, v (j)(v) 
e s ip'(u) + e s s jU Tp(u). 



Since R ;U v = R,vu (R( u ,v) E C 2 ), we obtain a first order PDE: 

(j)(v)s jV - ip(u)s )U + 4>'(v) - t//(u) = 0. 
The characteristic curves for eq. (117h are: 



v = —<t>{v) 

u = *p(u) 

s = 4/(v)-i//{u). 



(15) 



(16) 



(17) 



(18) 



We can find explicitly the functions <f>(v) and ip(u). In the standard gauge 
s{u,0) = so CES [[]) yields 



R, uv = F(R) = iP\u), 

t 

F(R(u', 0))du' + V(0). 



(19) 
(20) 



o 



Eqs. (|TT|) and (|T5j) in the standard gauge yield: 

" R du' _ f R F(R') 



i>{u)= F(R)—dR' + R, v (0,0) = - 

h(R(u,0)) - h(R a ) 



dR' + a 



+ a 



h(R a — bu) 



(21) 



(Here we use the definition of F(R) - eq. (|3|) and the value of h(R a ) is from 
eq. (ED). 
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(22) 



In the same way, on u = 

<f>(v) = [ R F(R)^-dR> + R, u (0, 0) = [ R ^ldR' - b = 

JRa dR JRa a 

h(R(0,v))-h(R a ) b _ h(R a + av) 
a a 

Now, we can rewrite the characteristic as 
'v = a h(R a + av) 

u = b^hiRa - bu) (23) 
s = h'(R a - bu) - h'(R a + av). 

Note the behavior of the function R(u,v) in the region (n > 0,v > 0): 
From eqs. f)15|) and (|22p we obtain an expression for Ry. 

R u = - a _1 e s /i(i? a + av). (24) 

We find out (from the form of the function h(R)) that R(u, v) is monotonously 
decreasing in regions I and II (Fig. [2]) on the rays v = const in consistency 
with the known expression for v = (eq. (jlip ). In an analogous manner, 
from eqs. f)15[> and (|2ip . we obtain an expression for R jV : 

R v = b' l e s h{R a - bu). (25) 

Thus, R(u, v) is monotonously increasing in region I on the rays u = const 
in consistency with the known expression for u = (eq. (jlip ). In region 
II R(u, v) is monotonously decreasing on u = const. We summarize the 
behavior of R(u,v) in Fig [2] below. 

Next we study the singularity of s. As we approach the point (2,1) in 
Fig. O we see that since h'(R 2 ) > 0, h'(Ri) < and h(R) > in region / 
(R(u,v) is monotonous), 

u, v — > + 

(the limit is from the positive side of zero). Thus, u and v attain their 
asymptotic values. Therefore from the last equation in (|23p we obtain: 

s -> h'{R 2 ) - h'(R x ) = const > (26) 

Then, s — > oo. However, as we approach the point (3,1) we have two generic 
cases: 

1. > k-i - s = const < while v — * + s — > — oo. 
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Figure 2: The sign of partial derivatives of R in the different regions of the flux- free 
solution in the standard gauge. The indications R\, R2 and R3 
label the corresponding characteristic lines. 

2. k\ > &3 - s = const > while v — > + =>■ s — > 00. 
From (|23p we may also consider the flow along the characteristics: 

du _ ah(R a - bu) 

dv ~ bh(R a + av)' ^ ' 

The direction of this flow is represented in Fig [3l 




Figure 3: The characteristic curves in the flux-free (linear) case in the standard 
gauge. 
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3 The general characteristic equations 



In this section we derive the characteristic equations in their most general 
form, i.e., for any perturbation. Furthermore, we show that the characteris- 
tic curves have physical meaning - they are the curves of constant R. Then 
we show that the parameter t of the characteristic curves also has a physical 
interpretation since it is invariant under gauge transformations. 



3.1 Derivation of the characteristic equations 

Similarly to (|17p . we can derive an equation for s(u,v) for any perturbation 
to the flux-free case. In this case, however, the equations are non-linear 
(and, also, non local) . 

Integration of (0 , [7]) gives: 



e- s R, u = I e- s( - u '' v ^(u')du' + (j)(v), (28) 
e~ s R v = e-< u ' v '^(v')dv' + ip(u). (29) 



The functions (f>, tp, ^ are determined from R(0, v) and R(u, 0) (in the 
standard gauge) as shown below. For convenience, we rewrite (128|29|) as 

Ru = e s fi{u,v), (30) 
R,v = e s h(u,v). (31) 

Differentiate (j30j) with respect to v, (f3T1) with respect to u, and equate 
R,uv = R,vu to obtain 

s, v fl ~ s,«/2 + fi, v ~ h,u = 0. (32) 
The equations of the characteristic curves for this equation are: 

'v = - fi{u,v) 

u = f2(u,v) (33) 

J = fl tV — f2,u- 

From daUlEII) and ([33]), we obtain: 

dv fi R tU - [ ' 

It follows 
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Proposition 3.1 A curve in the domain u>0v>0isa characteristic 
curve of \32\i if and only if R = const on the curve. 

As in the linear case, we can obtain from (|28M3ip : 

Ruv = e s s, v { re- s{u '' v) V(u')du' + cf>(v)) - e s ( P ' s. v e- s ^(u')du' - <f>'(v)) 
Jo Jo 

(35) 

Rvu = e s s, u ( f V e-< u ' v '^(v')dv' + (u)) - e s { P s, u e~ s ^(v')dv' - tf'(u)). 
Jo Jo 

(36) 

From (|35p on u = in the standard gauge (s(0,v) = 0) we obtain R uv = 
4>'{v). From ([36]) on v = in the standard gauge (s(u,0) = 0) we obtain 
R,vu = ip'( u )- Equating these expressions to eq. ([T]), we arrive at the same 
results for the derivatives of ip(u) and <f>(v) as in the linear case: 

- h'(R{u, 0)) = F(R(u, 0)) = rp'(u) (37) 

- h'(R(0, v)) = F(R(0, «)) = 4>'{v). (38) 

Note also that in the standard gauge the conserved fluxes (0 [7]) takes, 
on the axis u = 0, v = 0, the form: 

V(u)=R, uu {u,0) (39) 
= R m (0,v). (40) 



Combining the above results we can rewrite the characteristic equations (|33p 
explicitly (in terms of the data R(u,0) and R(0,v)), as 

= _ p e ~ s R, uu (u', 0)du> + fhf (R{0, v')) dv' 
u = fe~ s R >m (0, v')dv' - j u ti {R{u', 0)) du' 

J = h\R(u,6)) - h\R(0,v)) + P(e- s ) jV R, uu (u',0)du' - R tVV {0,i 
3.2 The physical interpretation of t 

The parameter t of the characteristic curves is related to the proper time in 
the same way as the coordinates r* or t are related to i10. The proper time 
r is defined by: 

dr 2 = —dudv. (42) 

R2 



1 t here is the parameter that we use in the usual black-hole metric, r* is the tortoise 
coordinate, defined by r* = f - dr lm for example in the Schwarzschild metric. 
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From eqs. ([33]) and (|30|) we obtain: 



u 



R „e 



(43) 



From (1311) and (1331) we obtain 



? R,u- 



(44) 



Combining the last three results, we find that the proper time r and the 
parameter t are connected through the scalar e~ s R u R v : 



2c 1^, 'U R ii 



R^ 



(45) 



Therefore the parameter t is a geometric quantity and not coordinate de- 
pendent. In the linear case we have: 



2h(R) 



R2 



(46) 



It appears to be like the connection between r and t on r = const, where 



t is the usual Schwarzschild coordinate (recall that 



h(R) 
2Ri 



2/ in the static 



solution), but here the characteristic curves are more general. For example, 
the relation (|46p is valid for curves of the type u = const and v = const even 
though it is not true for the Schwarzschild t coordinate. Such curves are the 
characteristics along the "horizons". Furthermore, on each characteristic 
curve R = const and therefore (f) 2 = const. Thus, the proper time is 
proportional to the parameter t, except for the horizons (roots of h(R)) 
along which f = 0). To conclude, in the linear case the parameter t is 
regarded as an extension of the Schwarzschild coordinate t to include the 
horizons as well. 



4 Discussion 

Let us summarize the features of the singularity obtained in the linear case 
{R,uu( u , 0) = R )VV (0, v) = 0). We have a singularity that develops from 
regular initial conditions. This singularity forms at finite u and v - a fi- 
nite time blowup. Its structure suits the structure of the RNDS spacetime 
(with an additional "horizon" because of the additional zero in h(R)). The 
singularity consists of a "point singularity" where s — > co and a ray of singu- 
larity emerging from it where s — > ±oo, depending on inequalities between 
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.V 



x Hi 



Hi 



Figure 4: The general singularity structure. The three horizons are displayed, 
denoted by Hi (i — 1, 2, 3). The singular inner horizon is displayed by a thick line 
emerging from the point singularity in the middle. 

the derivatives of h(R) at the horizons. The structure is illustrated in Fig. 
IH These inequalities and the resulting singularities of s correspond to the 
blue-shift\red-shift at the inner horizon. In addition, we found out that each 
horizon corresponds to a certain R value (which is constant on the horizon); 
In the linear case the horizons are characteristic curves of the first order 
PDE (|17[) according to Proposition 3.1. In the non-linear case the curves of 
constant R deviate from the horizons due to the perturbations. Therefore 
we will use in the non- linear case the notation Hi (i = 1,2,3) for the hori- 
zons. H\,H2, H% correspond in the RNDS model to the cosmological, event 
and inner horizons respectively. 

We reformulate below the main ingredient of the conjectural for the 
non-linear case (nonzero flux in both directions) : 

Conjecture 

For small enough non-linear perturbations the structure of the linear case 
is preserved. That is, a singularity point is formed from which the inner 
horizon emerges (H3 in Fig. [J). In addition, we expect the gauge invariant 
quantities to diverge at the inner horizon. Since they are related to the 
mass inflation phenomenon mentioned in the introduction, we expect the 
conditions for their divergence to be the same as in the mass inflation model 
of Brady and Poisson 

At the table below we summarize the divergence conditions of different 
quantities at the inner horizon based on [1] for the linear and non-linear 

2 The full form of the conjecture can be found in [I] 
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cases. 



The quantity 


Type of divergence 


Condition for 








the divergence 




linear 


non-linear 




s 


+00 


+00 


ki > k 3 




—00 


— CXD 


k 3 > ki 


e~ s R v 




+ OO 


fc 3 > 2k\ 


e R^ V R^ U 




+ OO 


k 3 > 2kx 



Table 1: The divergence of various quantities at the inner horizon H3 expected to 
be found in the linear and non-linear cases. 



h(R) as a saw-tooth function 

From now on we write the characteristic equations for a particular case of 
the function h(R). We choose h(R) to be a saw-tooth function as illustrated 
in Fig. The origin of the (u, v) axes is positioned on the curve R = R a . 
The region where R a < R(u,v) < R\ will be referred as region 1. The 
region where Rb < R(u, v) < R a is region 2. These, and regions 3 and 4 are 
displayed in Fig. [H 







k 2 s 














R, 


r\ 


► 



R 



Figure 5: The function h(R) as a u saw-tooth" function. 

For the saw-tooth function the evolution equation for s ([2]) turns to be 
in each region s U y — whose solution is 

s(u,v) = si(u) + s 2 (v), (47) 

where si(u) is an arbitrary function of u and 52(f) - a function of v. 
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Figure 6: The regions 1,2,3,4 in the solution - the forms of these regions are deter- 
mined by the form of the function h(R). The general directions of the characteristic 
curves in each region are indicated by arrows. 

5 Reduction to a system of ODE 

In this section we reduce our semi-linear system to a system of four first 
order ODE along a characteristic curve. First we begin the discussion by 
investigating the characteristic curve R(u,v) = R a . We choose the pa- 
rameter t along this characteristic curve to be the independent variable for 
the equations. Our semi-linear system (jl|2p takes the following form near 
R(u,v) = R a : 

j R,uv = e s [(h + k 2 )9{R - R a ) - k 2 ] 

\s,uv = e s (k 1 + k 2 )5(R- R a ), { SJ 

where 9(s) is the Heaviside step function (6(s) = 1 for s > and 9(s) = 
for s < 0). These equations are valid in the neighborhood of R = R a . In 
the standard gauge s = s\(u) is the solution in region 1 and s = s 2 (v) is the 
solution in region 2 - see flU}. | On the curve R — R a : 

si(u) = s 2 {v). (49) 

3 We assume here that there is no "hiding" of the singularity — namely, the curve Rb 
is monotonously decreasing (Fig. |5J . This assumption is included in the statement that s 
is only a function of v. In the next section we will formulate sufficient conditions for it. 
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Eq. (|48p takes the form R iUV = k\e Sl ^ in region 1. We integrate this 
equation with respect to u: 

pu 

R v = I kie^'Uu' + R v (0,v), (50) 



and then differentiate with respect to t. Using (|43I44|) we obtain the follow- 
ing equation for Ry. 

R v = hR, v -R, vv (0,v)e- s R, u . (51) 



We obtain an equation for R u by integrating eq. (08J) in region 2 with 
respect to v, and then differentiate with respect to t (using again (|43I44|) ) : 

R,u = k 2 R,u + R,uu( u i 

0)e~ s R )V . (52) 

We obtain the third equation by integrating eq. (|48|) with respect to u in 
region 2 (in the standard gauge s tV (0, v) = 0): 

s' 2 (v) = s >v = (k 1 + k 2 ) [ U e< u '^5(R-R a )du' = (h + k 2 ) , p ^ " > , , (53) 
where «o is the solution of R(uq,v) = R a . Then using (|44p we have 

'- (h + k -' ) W2kj\' (54) 

where we denote by R :U (R a ) the function R u on the curve R(u,v) = R a . 
Note that the equation for s is written in a form which is valid in all of region 
2. When we add ([43]) to the above three equations we get the complete 
system of ODE on the curve R = R a : 

R,v = k\ R jV — R jVV (0, f) e s R jU 
R, u = k 2 R, u + R,uu(u,0)e~ s R, v 

ii = e s R >v . 

Our goal now is to extend this set of four equations to any characteristic 
curve in region 2. It is only the first equation which has to be replaced in 
region 2. For this we recall the characteristic equations (|4ip and equate the 
equation for u along the characteristic with (|43H . After differentiation with 
respect to t and using (|4"3l |4*4"|) we obtain 

Rv = ~k 2 R, v - (h + k 2 ) R : u - R, vv (0, v) e~ s R, u . (56) 
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This set of four equations is accompanied by initial conditions on the u 
axis for (R iV , R tU , s,u) derived from our original semi-linear system. For a 
characteristic curve in the standard gauge originated at (u, 0) we obtain: 

(R jV (0)=R tV (0,0)-k 2 u 
R, u (0) = R >u (%0) 

a(0) = (57j 



Mo) 



u. 



In regions 3 and 4 we shall only need the equation for the function s. By 
the same argument as in the previous discussions, s(u,v) is only a function 
of v in region 4 (in the standard gauge). We denote s in region 4 by s±(v). In 
region 3 we can write the solution as s(u,v) = s±(v) + ss(u), where s^iv) is 
the solution in region 4 (by continuity) . Now we write our semi- linear system 
(Q3 [2]) with the relevant terms in the neighborhood of the curve R = R : 

R,uv = e s [-(k 2 + k 3 )6(R-R b ) + k 3 ] 

s,uv = e s [(h + k 2 )5{R - R a ) - (k 2 + k 3 )S(R - R b )}. {06> 



In regions 3 and 4 we have to integrate the second of (J58J) with respect 
to u over the delta functions across R = R a and R = R b and use the first 
equation in (|58l) to obtain: 

where R )U (R a ) is the function R tU (v) on the curve R = R a and R tU (Rb) is 
the function R u (v) on the curve R = R b . 



6 The non-linear case 

6.1 Singularity formation in the non-linear case 

The structure of Fig where R a is the boundary between regions 1 and 
2, implies i? n (0,0) < 0. Therefore there exists t > 0, so that i? n < on 
R = R a for < t < t. On this segment, (|55]) holds while the third equation 
is simplified into s = ki + k 2 . Then in the standard gauge we obtain: 

s = (k 1 + k 2 )t. (60) 

Proposition 6.1 Assume mi < i?„(0,v) < Mi and —M 2 < R :U (u,0) < 
—m 2 where mi, M\, m 2 , M 2 are arbitrary positive constants. Then, there 
exists singularity s — > oo on the curve R = R a at finite values of u and v. 
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Proof 

Using eqs. ((33} 133]) and the first two equations in ((33]) we obtain 
J2 tW (t) = fae* 1 * J*e- kltf R tV (0,v)dt? + R tV (0,v) 



(61) 



\i2,«(t) = he k2t jy k2t 'R, u (u, 0)dt' + i2, u (u, 0). 
From this and the assumption of the proposition on R v (0,v) and R u (u,0): 



mie klt < R )V (t) < Mie klt 
-M 2 e k2t < R, u {t) < -m 2 e k2t . 

Now we substitute OT) into (H3J HHD to obtain 



u = e-^ +k ^R tV 
v = - e -(ki+k2)t' R 



(62) 



(63) 



which, together with ((62]) . implies 

\^x#(i- e -n (64) 

In conclusion, when i — ► oo, s — > oo but lim^oo u(i) < 4p and lim^oo < 



&2 

Thus, the singularity forms at a finite u and v. □ 



Under the assumption of Proposition HTTl we obtain from (|62p that R jU (t) < 
for t G [0, oo) on R = R a . The same conclusion is valid for any character- 
istic curve in region 2 since the equation for takes the same form as on 
R = Ra- 
in Fig. [6] we assume that the curve R = Rb is monotonously decreasing. 
This is guaranteed under the following reasonable assumption 

Assumption 6.1 R w (0, v) < (the weak energy condition) 

The characteristic curve which plays an important role in determining 
whether Rb is monotonous or not is the curve which starts on the u axis 
at the point Uh = k^ 1 R )V (0,0). We denote the constant value of R on this 
curve by Rh- In the next result we show that R = Rh is monotonously 
decreasing, and that its Uh is the minimal value of u with this property. 

Proposition 6.2 Let Assumption \6. l\ The characteristic curve originated 
at Uh = k^ 1 R t v(0,fy on the u axis is monotonously decreasing (u < 0) 
provided Ub > Uh (see Fig. [6]). 
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Proof Let us write the characteristic equation for u(t) (|4ip in region 2: 

u= e' s R m (0,v')dv' + R V (0,0) - k 2 u. 
Jo 

Even though we still do not know the shape of the curve R = Rh , we 
know that in a small region near the axis u, s is a function only of v (a 
consequence of the choice in the standard gauge). This region is shown in 
Figured There exists < tf such that for any t £ [0, tf], s = s(v). Now we 
differentiate this equation with respect to t and use (|44p to obtain 

u + k 2 u + e ~ 2s(u) R, m (0, v)R >u = (65) 

By the choice Uh = k^ 1 R, v (0, 0) we have u(0) = 0. Next we show that 
u(t) < for any t £ (0,t/). Substituting u(t) = c(t)e~ k2t into §S§ we find 

c(t)e~ k2t = -e- 2s R >vv (0,v)R, u . (66) 

Since R >vv (0, v ) < (assumption l6.ll - the weak energy condition) and R tU < 
in region 2 (v = —e~ s R )U > in region 2 from (|44j) ). we find that c(t) < 0. 
In addition, c(0) = u(0) = 0. Therefore n < for any t E (0, t/). Note that 
v > in region 2 where s = s(u), thus ^ < on the characteristic curve 
R = Rh- 




Figure 7: The beginning of the curve R — Rh (the curve that starts at the point 
Uh = k^ 1 R, v (0, 0) ). Here we assume that the curve starts monotonously decreasing. 
According to Proposition \6.2\ this is the correct possibility. 

We now argue that tf = oo. Indeed, if tf < oo then, necessarily, u(tf) = 
0. This, however, is impossible via (|66p . whose right side is negative in 
region 2. □ 



18 




Figure 8: The curve R = Rh and some other characteristic curves. Notice that 
the curve R — Rh is the first characteristic curve which is monotonously decreasing 
(Erst in the sense of minimal u values) . The characteristic curves before R = Rh 
are not monotonous. 



The characteristic curve originated at (u/,,0) plays an important role in 
defining the boundary of region 2. If this point is part of the boundary of 
region 2, then the singularity is not hidden. Based on the assumption that 
— Mi < R tU (u,Q) < —m,2 we obtain: 

R a - M 2 u < R(u, 0) < R a - m 2 ii. 

Thus we can bound Uh from below: 

Ra — Rb 



u b > 



M 2 



(67) 



Therefore, a sufficient condition for u\, > u/ l (= k 2 1 R tV (0, 0)) in the non- 
linear case is 

R a - k 2 1 M 2 R, v (0,0) > R b . (68) 



Proposition 6.3 All the characteristic curves in region 2 meet at one point 
where s — > oo. 

Proof 

First we argue that any characteristic curve R ^ R a cannot intersect with 
R a at t < oo. If such an intersection takes place, then the gradient of R will 
blow up. The latter is impossible for t < oo according to ([62]) . The same 
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argument applies for any two characteristic curves in region 2. We conclude 
that the characteristic curves in region 2 cannot intersect for any t < oo. 

We showed in Proposition 16.11 that the singularity exists on the curve 
R = R a . Let us denote the point singularity on the curve R = R a by 
(u s ,v s ). Since in region 2, s(u,v) is a function only of v, we have the locus 
of the singularities in region 2: 

R°°(u ) = lim R(t) ; v(0) = s(0) = 0, it(0) = u , R(0) = R(u , 0), 

t— >oo 

where < no < ttfe (see Fig. [9] where the locus and one of the characteristic 
curves R = R are displayed). On R°°(u) we have a divergence of s via 
(|60p . Since no two characteristic curves can intersect at t < oo in region 
2, then R°°(ui) > R°°(u2) if < u\ < U2 < u^. If, for some tii < U2, 
R°°(ui) > R°°{u2), then 

f 2 R™du = i?°°(n 2 ) - R^iut) < , 

J Ul 

which implies that \R^\ < oo for almost any u in the interval (ui,^). But 
R^ = lim^oo R, u (t) = oo (from ([62]) ) where R, u (t) is given by (fBTj) . and we 
obtain a contradiction. It follows that is a constant in any such interval. 
□ 




Figure 9: The curve R = R and the curve R = R a hit the ray of singularity (R°°). 
The situation in this figure is impossible as we show in the proof of Proposition 8.3 

In order to proceed to regions 3 and 4, we need another result about 
region 2 which is formulated in the following Lemma. 
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Lemma 6.1 In region 2 n' U (R") ~^ 1 when t — ► oo, for any R' , R E [-R a ,i?b] 



Proo/ 

From equation (|48p we obtain, analogously to ([50 



R u = -k 2 /V^W + 22,„M) . 



This equation is valid for any characteristic curve in region 2. So, look at 
the limit 

iUi?) = -feJV^W + iU"'("),0) 
v-*vsR >u (R») v^v s -k 2 fe s *( v ')dv' + R }U (u"(v),0)' K ' 

u'(y) and u"(v) describe the curves R(u,v) = R! and R(u,v) = R" respec- 
tively. According to (|60p the integrals in the denominator and the numerator 
diverge. Using L'Hospital's rule we obtain: 

lim **M = hm -^Wl-^"^ (70 ) 



]im v -+ Vs e S2 ^j^ = on any characteristic because lim v ^ Vg J v e S2 ^ v '^dv' di- 

v ° R,u(R") 



verges but lim^^ f V ^jr^ v ' = u ( v s) is finite. Therefore ]hn v - tVa ^' u mJ\ = 1- 



□ 

Let us summarize the properties of the inner horizon of our toy-model 
in the following proposition. 



Proposition 6.4 Assume t68\) . Then, at the inner horizon H%: 

a. If k^ > ki, s — » oo, while if ki < k^, s — > — oo. 

b. R )V e~ s diverges if k% > 2k\ and converges if k% < 2ki . 

c. e~ s R )V R^ u diverges if k$ > 2k\ and converges if k^ < 2k\. 

Proof 

a. In region 4 we have an expression for s±(v) as a function of s 2 (v) and 
R :U in region 2 (eq. (|59|) ): 



s 4 (u) 



1 



R,u(Rb 



k2 + h _ (h + k 2 ) R, u (Rb) 
R u [R a ) 



s^'W. (71) 



Let us look at the integral when v — > v s . The integrand is written in 
terms of quantities in region 2. Recall that in region 2 R tU < 0. Then since 
s 2 (v) — > oo when v — > v s and using Lemma 6.1 we see that if ks > k\ then 
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S4,(v) — > — oo and if k% > k%, s^(v) — > 00. Note that S4(i>) is regular in all 
region 4 for v < v s . 

In region 3, although s(u, v) = 54(7;) + S3(u), we have the same type 
of singularity as in region 4 along the horizon H3 (except for the point 
singularity). S3(u) is regular in region 3, except for the point singularity. 
The reason for this regularity of S3 (u) is the following: If it was not regular in 
region 3, we would have a singularity in the same u on the curve R = Rb (Fig. 
[TO]) . This, of course, contradicts our previous results about the regularity 
of s in region 2. Therefore the only point where 53(11) might diverge is the 
point singularity. 




Figure 10: A situation where there is a singularity in S3(u) in u s > u s . Then we 
would have a singularity on R = R a at a regularity point of S3 . 

b. From the inequality for R u (|62p . which is valid in all region 2, we 
obtain the asymptotic behavior i? n ~ e h2t . Substituting the asymptotic 
expressions of R^ u and e S2 ^ into eq. (|44|) . we find that v ~ e~ klt . Let us 
denote by Av the distance to the singularity on the v axes. Then, asymptoti- 
cally, Av ~ e _fcl * near the singularity (t — > 00). Substituting the asymptotic 
expressions of R :U , e S2 ^ and v into eq. (HU we obtain: 

Substituting Av, we obtain: 

e S4 ~ Aw^t" 1 (73) 
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in region 4. Equating the expression for ii in (|4ip with (|43p and substituting 
the asymptotic behavior (|73p we obtain 

(Av') 1 '^ R vv {0,v')dv' - h'{R(u',0))du'. (74) 

Therefore e~ s R )V diverges near the inner horizon H% if k% > 2k\ and con- 
verges if £?3 < 2fei. 

c. Recall the evolution equation for R(u, v) in region 4: 

R,uv — 4 ^ 

Integrating this equation with respect to v we obtain: 

R, u = h J e s ^'W + R )U (u, 0). (75) 

Substituting the asymptotic expression for e Si (|73p into the expression for 
R u , we find that at the inner horizon R u converges to a constant different 
than zero. Therefore the divergences of e~ s R )V R )U and e~ s R tV at the inner 
horizon are the same. □ 



6.2 The No-Hair Phenomenon 

Our semi-linear system demonstrate a no-hair phenomenon. Asymptotically 
near the singularity the contribution of the perturbations disappear. The 
best way to see it is to look at the equations and the asymptotic form of their 
components expressed with the parameter t. This parameter is appropriate 
for this purpose as well since t — > oo when we expect the perturbations 
to vanish. Considering (|55p and noting that the equation for s can be 
approximated by s ~ k\ + k^, we obtain the asymptotic forms of e s and 
R )U in region 2 as R yU ~ e k2t , e s ~ e (k 1 +k 2 )t^ Substituting these asymptotic 
forms in the perturbative term in the first equation of (|55p we obtain 

R, vv (0,v)e- s R, u ~ R, vv (0,v)e- klt . 

Thus, when t — > oo, the contribution of the perturbation decreases expo- 
nentially. Thus the asymptotic form of the equation for R tV is R )V ~ k\R )V . 
Therefore the asymptotic form of R^ v is R^ v ~ e klt . Substituting the latter 
in the second equation of (|55p we find that: 

R, uu (u, 0)R, v e~ s ~ R, uu (u,0)e- k2t . 

Hence, when t — > oo, the effect of the perturbations on the equations (and 
their solution) vanishes. This is the no-hair phenomenon - the characters of 
the black hole do not depend on the explicit form of the perturbations. 
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7 Conclusions 



The main goal of this research was to provide evidence for the validity of the 
conjecture, which was described in section 0] (and in [lj). We actually proved 
the conjecture for the particular case of h(R) as a saw-tooth function under 
few assumptions. The saw-tooth function is assumed to be generic enough, 
so any other function h{R) will yield a similar structure of singularity for 
an open neighborhood of the linear case in the space of initial data. 

For h(R) a saw-tooth function our semi-linear system demonstrates a 
finite time blowup in region 2 under the following assumptions: 

• mi < R <v (0, v) < Mi and — Mi < R lU (u, 0) < — mi where mi, Mi, mi, Mi 
are arbitrary positive constants. 

• R <vv (0,v) < (the weak energy condition). 

. Ra _M^m >Rb . 

In addition our semi-linear system demonstrates a no-hair phenomenon in 
region 2. We proved that the singularity in region 2 is actually a point 
singularity (which corresponds to future timelike infinity). The asymptotic 
expressions near the singular lines in regions 3 and 4 are in full agreement 
with the conjectured asymptotic behavior (Table 5.1): 

M_i 

• e s ~ (Au) fc i 

• e' s R, v R iU ~ (Av) 2 ~^ . 

The "half gauge invariant scalar" e~ s R >v is related to the divergent flux 
at the inner horizon in various models (see [1]). The gauge invariant scalar 
e~ s R jV R^ u is related to the mass function and the scalar curvature singularity 
which is formed when the model contains two intersecting streams of null 
fluid. 

The curve R = Rh (the curve that starts at the point Uh = k^ 1 R, v (®, 0) 
on the u axis) has turned out to be very important in the definition of region 
2 since it is the first curve in the region which is monotonously decreasing 
(see Fig. [HJ). This curve has an important physical meaning for our model: 
It is the first characteristic curve (R = const) which is completely inside 
the trapped region. Namely, R tV ,R )U < along this curve. The apparent 
horizon is defined as the outer boundary of the trapped region inside the 
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black hole. In the linear case (which corresponds to the structure of a static 
black hole) the event horizon and the apparent horizon coincide (see Fig. 
[2]). In the non-linear case the apparent horizon is inside the black hole. In 
Fig. [8] we can see how the apparent horizon passes if we join the extremal 
points of the curves behind the horizon H2 by a line - the parts of the curves 
which are monotonously decreasing are inside the trapped region. 

In conclusion, we found that our semi-linear system with h(R) as a saw- 
tooth function demonstrates many features that have physical meaning and 
are analogous to the black-hole inner structure (the causal structure and 
location of the singularity, the infinite blue shift and the inner horizon, mass 
inflation, the trapped regions etc.). We even found in the model features 
of the spacetime outside the black-hole (e.g. no-hair principle). Further 
research is needed to confirm (or contradict) the conjecture under general 
data. Even more challenging question is the validity of the conjecture for a 
general h(R). 
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